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PREFACE 


In  1973,  the  Aerospace  Applications  Studies  Committee  (AASC)  of  the  Advisory 
Group  for  Aerospace  Research  and  Development  (AGARD),  North  Atlantic  Treaty 
Organization  (NATO),  sponsored  a study  on  the  application  of  night  vision  devices 
to  fast  combat  aircraft.  During  the  study  it  became  apparent  that  the  assumed 
weather  conditions  — highly  smoothed  10-year  averages  — were  far  too  uniform  to 
give  realistic  results.  Curiosity  about  the  variations  of  unsmoothed  weather  data 
led  to  a proposal  to  the  AASC  by  L.  M.  Biberman  of  the  Institute  for  Defense 
Analyses  (IDA)  and  M.  H.  A.  Deller  of  the  Royal  Aircraft  Establishment  (RAE), 
Famborough,  that  the  problem  be  investigated  in  some  detail  to  learn  the  effects 
of  terrain  masking,  cloud  obscuration,  and  hour-by-hour  weather  variations  at  a 
number  of  European  locations. 

The  resulting  study,  ^ published  in  five  parts,  contains  estimates  of  the  hourly, 
daily,  and  seasonal  effects  of  the  actual  weather  at  Hannover,  Federal  Republic  of 
Germany,  in  1970  on  the  performance  of  electrooptic2d  imaging  sensors.  The  questions 
we  hope  to  answer  are  how  great  these  effects  tire  and  when  and  how  often  they  occur. 

Part  1 of  the  study  discusses  methodology  and  samples  the  results  of  calculations. 
Part  2,  in  another,  classified  volume  (IDA  Paper  P-1124),  presents  complete  results 
for  a forward-looking  infrared  (FLIR)  device  in  the  8.  5-11  /um  bemd  and  analyzes  the 
impact  of  weather  on  operations  and  operational  planning.  Part  3 (IDA  Paper  P-1128) 
compares  FLIR  performance  in  the  3.4-4.  2 Mm  and  8.  5-11  Mm  bands.  Part  4 (IDA 

^"Effect  of  Weather  at  Hannover,  Federal  Republic  of  Germany,  on  Performance  of 
Electrooptical  Imaging  Systems",  References  1-5. 


Paper  P-1202)  reports  on  calculations  for  active  and  passive  television.  P^u:t  5 (IDA 
Paper  P-1203)  compares  the  performance  of  active  television  and  several  different 
FLIRs. 


This  note  was  written  as  an  explanatory  addendum  to  the  Study  and  thus  it  bears 
the  same  title.  The  note  documents  the  computer  program  (called  Program  FLIR) 
to  calculate  the  probabilities  of  detection  and  recognition  of  a target  by  an  observer 
using  a FLIR  sensor.  It  was  written  to  summarize  the  basic  concepts  behind  the 
calculation  procedures  in  Program  FUR  and  to  outline  those  procedures.  For  more 
details  about  the  physics,  bar-pattern  criteria  and  role  that  weather  plays  in  the 
calculations  refer  to  Parts  1-5  of  the  Study  mentioned  above. 
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I.  SUMMARY  OF  PROGRAM  FLIR 


Program  FLIR  was  written  to  calculate  the  probabilities  of  detection  and  recogni- 
tion of  a given  sized  target  by  an  observer  employing  an  8.  5-11  Mm  forward-looking 
infrared  (FLIR)  sensor.  The  probabilities  are  calculated  for  various  ranges  given 
transmission  at  those  particular  ranges.  A basic  outline  of  the  program  follows 
including  necessary  inputs,  types  of  calculations  performed  and  the  output. 

The  input  of  Program  FLIR  includes  specifying  the  base  field  of  view  (FOV)  of  the 
sensor,  scaling  factors  for  subsequent  FOV's,  and  the  appropriate  curves  which  re- 
late spatial  frequency  to  minimum  resolvable  temperature  (MRT)  and  characterize 
the  design  of  the  sensor  (Figure  1).  Target  characteristics  to  be  input  are  critical 
dimension,  differential  temperature  ( AT)  between  the  target  and  its  background  and 
aspect  factors.  The  ranges  for  which  probabilities  are  to  be  calculated  are  also 
entered.  Most  of  the  data  input  to  program  FLIR  consists  of  transmission  values 
which  were  previously  calculated  and  written  onto  tape.  The  transmissions  must  be 
computed  for  the  ranges  specified  as  input  to  Program  FLIR  and  are  usually  calculated 
hour  by  hour  for  a particular  month,^  The  transmissions  are  based  on  the  weather 
(such  as  air  temperature,  dew  point  and  visibility)  recorded  for  those  hours  for  that 
month  for  the  location  under  consideration.  One  method  for  obtaining  the  transmission 
from  weather  statistics  is  by  using  the  model  Lowtran  3B,  an  atmospheric  trans- 
mittance model  developed  by  the  Air  Force  Geophysics  Laboratory  (Ref.  6). 

^Lowtran  3B  includes  several  additions  and  updates  to  the  Lowtran  3 and  3A  models  re- 
ported in  Ref.  1.  The  major  additions  are  the  inclusion  of  water  vapor  continuum 
attenuation  in  the  3.  5 to  4.  2 Mm  region,  and  a temperature  dependence  to  the  H2O  con- 
tinuum attenuation  coefficient  in  both  the  4 Mm  and  10  Mm  regions. 

2This  data  becomes  columns  1 and  2 in  Table  2. 
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The  calculations  performed  by  Program  FLIP  using  the  input  variables  are 
summarized  here.  For  greater  detail  refer  to  Chapter  III,  Users  Guide  to  Program 
FLIP  and  Chapter  IV,  A Listing  (with  documentation)  of  Program  FLIP. 

First  a subroutine  (Spatial)  calculates  the  angular  subtense  of  the  target  (in 
milliradians)  and  returns  the  values  of  the  spatial  frequencies  in  cycles/milliradian 
for  both  detection  and  recognition  at  a given  range.  Given  the  computed  spatial  fre- 
quency for  each  range  another  subroutine  (Pesolve)  scales  it  according  to  desired  FOV. 
The  scaled  spatial  frequency  is  then  used  to  find  the  corresponding  MPT  value  by 
linear  interpolation  between  points  on  the  MPT  curve  input  to  the  subroutine.  This 
approximation  is  used  when  the  MPT  value  is  not  available  for  the  FOV  of  interest, 
but  is  available  for  an  FOV  of  relatively  close  size.  The  MPT  value  determined  is 
then  corrected  for  aspect  ratio  of  the  target  (Pef.  1,  Appendix  D,  "First  Order 
Corrections  to  Bar-Pattern  Data"). 


Next,  by  hour  and  by  r<mge  the  transmission  data  is  read  in  and  the  apparent  target 
temperature  ( ATgpp)  is  Ccilculated  by  multiplying  the  transmittance  ( ’’atm) 
of  the  target  which  has  been  input.  ^ Finally,  given  the  ratio  AT^pp/MPT  corrected 
for  aspect,  which  is  the  normalized  displayed-signal-to-noise  ratio,  the  probability 
of  detection  (or  recognition)  is  determined  by  calling  the  last  subroutine  (Cuprob). 

This  subroutine  determines  the  probability  of  looking  up  the  signal-to-noise  ratio  on 
a cumulative  normal  probability  curve  where  the  value  ATjpp/MPT  = 1.  0 corresponds 
to  a probability  of  detection  of  50  percent. 


= AT  • 


^atm* 


The  output  of  Program  FLIR  is  a computer  tape  which  has  written  on  it  tables 
which  look  like  Tables  1 cuid  2.  Table  1 is  a header  table  to  the  FLIR  tape  and  con- 
tains the  MRT  data  for  each  range  and  FOV.  Table  2 contains  the  probability  of  de- 
tection and  recognition  data  for  two  hours  by  range  and  FOV.  There  is  a table  for 
every  hour  of  a month  on  a typical  output  tape.  Using  this  FLIR  output  tape  various 
plots  and  other  outputs  can  be  made  (see  Chapter  V.  Examples  of  Outputs  Which  Can 
Be  Generated  From  the  Output  Tape  of  Progrcun  FLIR). 
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OUTPUT  FROM  PROGRAM  FLi 


Table  2.  EXAMPLE  OF  HOURLY  TABLES  OUTPUT  BY  PROGRAM  FLIR 


r 

IL  CAVEATS 

' This  section  reiterates  the  caveats  regcurding  the  data  produced  by  Program  FLIR, 

i 

the  input  data  which  had  been  calculated  by  Program  LOWTRAN  and  the  masking  effects 
as  discussed  in  Part  1 of  this  study. 

I 1.  Data  Calculated  by  Program  FLIR 

i 

The  data  computed  by  this  program  are  necessary  but  insufficient  for  assessing 

I 

I the  overall  effectiveness  of  electrooptical  imaging  sensors  aboard  attack  aircraft 

1 

used  against  ground  targets.  It  must  be  recogmzed,  however,  that  the  data  will  almost 
^ always  represent  an  upper  bound  of  performance,  since  the  computed  probabilities  of 

detection  and  recognition  at  given  ranges  assume  that  the  observer  is  already  looking 
at  the  area  of  the  display  that  coincides  with  the  position  of  the  target  within  the 

! 

field  of  view.  Realistic  assessments  of  the  observer's  capabilities  will  require  better 

1 data  than  are  currently  available  on  his  display  search  time,  on  his  dynamic  task 

i 

I 

i performance  (including  his  target  recognition  and  weapon  aiming  time),  and  on  the  de- 

; gradations  to  be  expected  from  both  the  airborne  environment  and  the  actual  opera- 
tional environment.  ^ 

: < 

■ 

No  matter  how  good  the  viewing  conditions  and  equipment,  a drowsy  or  disinterested 
observer  will  not  do  very  well.  It  was  not  within  the  province  or  competence  of  this 
study  to  ascertain  motivation  or  interest  on  the  part  of  the  observer.  Our  calculations 
are  based  on  more  than  200, 000  data  points  for  the  performance  of  serious  observers 
looking  at  a variety  of  targets  displayed  against  various  noisy  backgrounds  but  we  do  not 
know  how  much  to  degrade  our  results  to  cover  various  tactical  situations  or  observers 
who  are  not  very  attentive. 
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We  have  <dso  excluded  degradations  due  to  exposure  of  airborne  observers  to 
buffeting  or  g-loading.  Preliminary  experiments  completed  in  M«y  1975,  ^ showed 
that  both  buf feting  and  g-loading  degrade  observer  performance,  but  no  analysis  of 
the  frequency  and  severity  of  these  effects  suitable  for  use  in  modifying  our  results 
is  yet  available. 

Arguments  about  the  modeling  of  recognition  range  and  about  how  to  define  "identifi- 
cation" remain  unresolved.  Semantics  gain  in  importance  as  the  tactical  problem  shifts 
from  detection  to  identification.  In  military  operations,  recognition  of  an  electro- 
optical  image  of  a target  is  very  closely  related  to  circumstances.  Given  the  appropri- 
ate background  intelligence,  a senor  operator  can  positively  translate  a series  of 
unresolved  specks  moving  along  a road  into  a column  of  trucks  or  tanks,  poor 
optical  quality  notwithstanding. 

For  detection  of  tactical  vehicular  tugets,  we  use  the  criterion  of  two  lines^ 
across  the  minor  dimension  of  the  tcurget,  on  which  there  is  generally  good  agreement. 
For  recognition  we  use  four  resolvable  line  pairs  across  the  minor  dimension  of  the 
target.  In  undemanding  situations,  some  people  elect  to  use  three  line  pairs  as  a 
criterion  for  recog-dtion.  In  bad  clutter,  some  use  a criterion  of  four  and  a half  or 
five  line  pairs,  but  we  prefer  to  stay  with  four  line  pairs  and  to  increase  the  signal- 
to-noise  threshold  (Refs.  7,  8). 

In  the  computational  program  presented  herein  we  treat  in  depth  only  data  on 
FLIR. 

^Discussed  in  Appendix  H,  Reference  1. 

^Two  lines,  one  line  pair,  and  one  cycle  are  synonymous. 
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2.  Computation  Model  Lowtran 


In  our  computations  we  were  bothered  by  four  problems: 

1.  The  overly  pessimistic  predictions  of  Lowtran  3 for  the  water  vapor  continuum. 
We  have  solved  this,  bringing  calculated  results  into  line  with  measured  values, 
by  correcting  the  Lowtran  predictions  for  the  continuum  and  its  temperature 
coefficient.  The  Lowtran  36  version  incorporating  this  change  is  now  being  used 
(Ref.  6). 

2.  The  weakest  part  of  Program  FLIR  is  the  aerosol  model  used  in  Lowtran  for  l,ow 
visibility  conditions.  Existing  Lowtran  aerosol  models  can  yield  a large  variation 
in  computed  detection  and  recognition  range  for  a given  FLIR.  Our  best  recom- 
mendation for  now  is  to  use  the  Lowtran  36  maritime  aerosol  model  for  Central 
Europe.  We  at  IDA  are  attempting  to  develop  suitable  subroutines  for  aerosols 
elsewhere  in  the  world  and  will  publish  these  as  they  become  available. 

3.  The  vertical  lapse  rate  for  mists  and  fogs.  The  craiditions  on  which  we  have  based 
our  calculations  are  valid  for  ground-to-ground  observations.  If  an  airborne 
sensor  is  looking  down  from  200  feet  above  ground,  however,  and  if  there  is  a 

fog  l^er  100  ft  thick,  the  path  through  fog  is  only  half  what  we  have  used  in  our 
computations.  At  present,  we  have  almost  no  data  on  the  layering  of  fog  and 
haze  at  Hannover  or  anywhere  else. 

4.  6y  international  convention,  visibility  exceeding  10  km  is  reported  as  infinite 
in  aviation  weather  data.  We  have  examined  the  effect  of  truncating  visibility 
at  10  km  by  recomputing  for  20  km.  Almost  no  change  in  the  statistics  could  be 
found,  since  visibility  dominates  only  when  its  values  are  small. 

3.  Masking  Effects 

The  effects  of  cloud  masking  and  terrain  masking  are  not  included  in  our  models 
of  probability  of  detection  or  recognition  but  must  be  considered  in  operations. 
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m.  USERS  GUIDE  TO  PROGRAM  FLIR 


Program  FLIR  is  written  in  FORTRAN  designed  to  r\in  on  a CDC  6400  computer. 
The  deck  is  punched  BCD  (on  a 026  IBM  keypunch).  Tape  1 u*  read  as  input,  Tape  2 
is  written  as  output.  Line  numbers  referred  to  in  the  following  sections  correspond 
to  numbers  assigned  to  each  line  of  Program  FLIR  and  its  subroutines.  The  lines 
and  numbers  are  listed  in  the  next  section  of  this  report.  Chapter  IV. 

A.  INPUTS  AND  FORMATS 
1.  Inputs  by  Data  Card 


Format 

Variable  Name 

Description  of  Variable 

Card  1 
(Line  41): 

IS 

NMO 

No.  months  for  which  FLIR 
is  to  be  run. 

F5.0 

FOVl 

Base  field  of  view. 

IS 

NF 

No.  mult^les  of  FOVl  to  be 
considered. 

FS.  0 

RMAX 

Maximum  range  in  km  - can  be 
either  10.  or  20.  (If  RMAX=10. 
range=.  S to  10.  km  in  steps  of 
. S km.  If  RMAX=20.  range=l. 
to  20.  km  in  steps  of  1 km). 

Card  2 
(Line  61): 

11 

ICURVE 

MRT  curve  selector;  MRT  values 
are  input  in  data  statements  in 
Subroutine  Resolve;  ICURVE -1 
unless  more  than  one  MRT  curve 
is  listed  in  Resolve  - if  more  than 
one,  ICURVE  equals  the  no.  the 
MRT  curve  desired  occupies  in 
the  list  of  MRT  curves  (see  next 
Section). 

4X 

FIO.O 

SIZE 

Minor  dimension  of  target  (in 
meters). 

FIO.O 

TARGT 

Temperature  (deg.  K)  diEf  eren- 
tial  of  target  from  background 
(AT). 
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Format 


Variable  Name 


FIO.O 

FIO.O 


ASPECT  (1) 
ASPECT(2) 


Description  of  Variable 


£/7  for  detection. 


^t/7  for  recognition.  I is  as- 
pect ratio  of  one  bar  in  the 
equivalent  bar  pattern  (Figure 
28,  p.  51,  Ref.  1).  For  front 
aspect  tank  detection  the  tank 
is  about  square  and  the  bar  as- 
pect is  due  to  one-half  a square 
or  2:1  aspect.  For  recognition 
there  are  four  line  pairs  so  one 
bar  represents  1/8  of  a square 
or  8:1  aspect  ratio. 


2.  Input  by  Data  Statement 

The  MRT  curve  points  are  input  by  a data  statement  in  tlie  program  deck  in  Subroutine 
Resolve  (lines  27R  - 36  R).  The  following  is  an  example  of  the  data  statement  when  it 


contains  five  curves  at  the  same  time. 


para  Cmpvf/1  . t2,  »3.  »4.  *7*  »8*0»  » T-1  X 

• I Y 

1 • • 2» 1 3« • 4 • » 5« * 6 • * 7 • • • 7*0 • * T*2  X 

.nO7.,dl8.,03**i06t.lOZ»,Z«.*.l.*7*0.*  T^'v 

1 • » 2»  » 3«  * ♦ « * 5*  ♦ 6 • 1 7 • • • 9*  • 6*<' « * T«3  x 

.nO2,THO42».00B-».Oi3..D26f .Oe ?77fffl77r7irr;s2'.'^ *TTri  T-?  Y 

2«»4»t6,»«,»10»»l0»0»» H«1  X 

•0062t»014**036t*05^*»0®H*1fl*rt,*  H*1  Y 

1 • > 1 • 8 1 2 • • 2 • 4 * 2 • 8 * 3 • 0 1 ^ • 2 * 3 • a « 3 • 6 • 3 • 8 1 3 • 9 1 3 • 99  ( 3*0  • t POO  x 

' »o6» • 16t  »3t  *42* ,b*  abt • • B4  fl •0^*l*24fl«32«3*o»/  PPD  Y 


Five  is  the  maximum  number  of  MRT  curves  that  may  be  entered  at  one  time. 
However,  only  one  curve  is  used  per  run  of  Program  FLIR.  The  curve  to  be  used  is 
specified  by  the  variable  ICURVE  (see  Input  Card  1 above).  In  the  excunple  shown 
abo%8  ICURVE  =1  would  select  MRT  curve  T-1,  ICURVE=5  would  select  MRT  curve 
PDD,  where  T-1  and  PDD  are  our  code  names  for  actual  FLIR  equipments. 
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There  are  two  lines  of  data  punched  per  MRT  curve  entered.  The  first  line 
labeled  "X"  contains  the  X -coordinates  of  the  curve  which  are  the  spati£d  frequency 
values  in  cycles/mrad.  The  second  line  labeled  "Y'*  contains  the  Y-coordiiuites  or 
the  MRT  values  (see  Figure  1 in  Summary). 
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3.  Data  Input  by  Magnetic  Tape 


The  previously  calculated  transmissions  are  read  from  tape  (Tape  1 on  program 
Ccurd).  The  transmission  tape  contains  one  file  per  month  to  be  processed.  Usually 
a maximum  of  two  files  are  on  one  transmission  tape.  Each  month's  file  contains 
three  lines  of  heading  or  title  information  to  be  skipped  over  when  being  read. 

Format  105  in  the  program  does  the  skipping: 

Line 

Read(l,  105)IBLANK  132 

105  Format  (//AlO)  133 

The  headings  are  followed  by  one  line  for  each  hour  in  the  month.  Each  hour's 
record  consists  of  3 integer  values  representing  month,  day,  and  hour,  and  20  real 
values  representing  the  fractional  transmittance  for  ranges  of  . 5 to  10.  km  in  steps 
of  . 5 km  or  1.  to  20.  km  in  steps  of  1.  km.  RMAX  on  input  card  1 above  indicates 
which  ranges  are  on  the  tape.  The  read  statement  for  one  hour' s data  in  this  form  is 

Line 

Read  (1, 1000)MDH,  TRANS  154 
1000  Format  (212, 14, 2X,  20F  6.3)  155 

If  a transmission  tape  of  another  format  is  used  the  two  read  statements  and 
formats  mentioned  here  will  have  to  be  revised.  If  different  ranges  are  to  be  used 
and  the  number  of  ranges  are  changed,  dimension  statements  may  have  to  be  changed 
along  with  DO  LOOP  counters  (e.  g. , Line  80  is  now  set  for  20  ranges)  and  scaling  of 
range  values  (e.  g. , Line  81). 
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B.  PROGRAM  FUR  DIVroED  INTO  SEGMENTS 


Lines 
34-66 
70-74 
r 80-123 


Input  cmd  printing  o£  2 data  c£u:ds. 

Write  headings  on  output  tape  (Tape  2). 

Do  Loop  5 is  executed  once  for  each  range. 

81  Define  range  for  this  execution  of  loop. 

87  Call  Subroutine  Spatial  which  calculates  angular  subtense 

of  target  (in  mrads)  and  the  spatial  frequency  (cycles/ 
mrad)  of  the  detection  and  recognition  criteria. 


“ 95-116  Do  Loop  2 is  executed  once  for  each  FOV. 

plOO-116  Do  Loop  2 is  executed  once  for  detection  and  once 
for  recognition  for  each  FOV. 

106  Call  Subroutine  Resolve 

27R-36R  Input  MRT  ciurve  X and  Y 

coordinates  in  data  statement. 

44R  Scale  spatial  frequency  calculated 

in  Subroutine  Spatial  for  the  FOV 
being  considered  in  this  pass 
through  Loop  2. 

S1R-83R  Linear  interpolation  along  MRT 
curve  to  find  correct  MRT  value 
corresponding  to  the  scaled 
spatial  frequency. 

112  Correct  MRT  value  returned  from  Subroutine 

Resolve  for  aspect  ratio. 


121  Write  table  as  heading  to  output  tape.  One  line  per  range  con- 

tains the  following  data:  angular  subtense  of  target,  spatial 
frequency  (cycles/ mrad)  for  detection,  MRT  values  for  detec- 
tion and  recognition  for  the  fields  of  view  considered. 
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I—  131-213  Do  Loop  50  is  executed  once  for  each  file  (or  month)  to  be  written 
on  Tape  2. 

p 154  Read  transmission  tape,  one  hour's  data  at  a time. 

163-170  Write  headings  on  Tape  2 for  each  hour's  table. 

174-209  Do  Loop  30  is  executed  once  for  each  r<mge. 

175  Define  the  range  for  this  execution  of  loop. 

179  Ccdculate  apparent  target  temperature. 

p 183 -205  Do  Loop  20  is  executed  once  for  column  of  pro- 
babHitiei  (number  of  columns  equals  2 times 
number  of  FOVs). 

196  Calculate  signal -to -noise  ratio. 

201  Call  Subroutine  Cuprob 

15C  Input  cumulative  normal 

probability  curve  coordinates. 

34C-37C  Interpolate  between  points  on 
noimial  curve  to  determine 
probability  of  detection  (or 
recognition)  corresponding  to 
calculated  signal -to -noise 
ratio. 

207  Write  on  Tape  2 one  line  per  rcmge  with  the  follow- 

ing information:  range,  apparent  target  tempera- 
ture, (probability  of  detection,  probability  of 
recognition  for  each  FOV). 

^ 211  Go  back  to  154  to  read  transmission  for  next  hour. 

214  STOP 


- 18  - 


»r'0»rir»rir>r»iririr»r  o tnnrtr  nricn  tr.trirn  trtcjrirn'.trjr.trtr.trir  eurtcirtr'rmrir.trtrtc  n 


IV.  LISTING  OF  PROGRAM  FLU? 


P90rt»*M  PL T R < INPUT t OUTPUT t TAPt 1 • T*PF2 ) 

program  FLTR  has  HRITTtN  AS  A BASIC  MOOCL  UP  THE  PORWaRO.LOOKINQ 
TNPRaREO  (PkIR)  DEVICE,  the  DEVELOPMENT  OP  THE  MODEL  AND  RESULTS 
OrTAIHEO  prom  its  use  are  RCPORTEO  IN  IDA  paper  P«Il?3t  EFFECT  OF 
HEATHER  AT  HANNOVER.  PMfl,  ON  PERPOhmANCE  OP  ELECTROOPTICaL  IMAGING 
systems.  PaHT  1 (REPERtO  TO  IN  THE  DOCUMENTATION  OF  THIS  PROGRAM 
A$  arFrOrt*) , 

program  FLIR  has  HRITTEN  to  read  an  HOUR  ST  HOUR  tape  OF  TRANSMISSIONS! 
TO  SCALE  MRT  values  FOR  DETECTION  «NO  mECOGNITIOn  AT  SELECTED 

Ranges  for  selected  fields  of  vie*  <fov»  ano  to  hrite  them  in'table 
Form  at  the  reginning  op  the  plir  taRE»  to  calculate  hour  bt  hour  for 
EACH  Range  and  For  each  fqv  the  apparent  target  temrfhature  ano 

PROBABTLITV  OF  DETECTION  ANO  RECOGNITION.  T» IS  PROGMAM  MAT  BE 

spgmenteo  into  three  Parts— 

PART  1.  ANRUT  VALUES  NECESSARY  PoR  EXECUTION  OF  THE  PROGRAM. 

PART  H.  LOOP  S— calculates  THE  MMT  values  and  WRITES  THEM  AS  A 
header  table  TO  THE  FLIh  TAPE. 

part  S.  loop  5o»CALCULATES  ano  mRITES  UN  THE  flip  tape  the  apparent 
TAHoET  temperature  ano  rrosarilities  one  pile 
(A  MONTHS  data  usually  EQUALS  ONE  FILE)  AT  A TIME 
HOUR  BY  hour. 

FOR  MORE  details  see  LINE  ST  LINE  UOCUMENTATION  OF  The  PROGRAM. 

OIMfmSIOM  ISTOP(IO) .FOV(Sl .ASPECT(2) .CPM(2) .RTMP(i0.20) 
dimension  mOH(3) .TRanS(20) .PRoBIIQ) 
real  Mrt 

TSTOP  hill  contain  PILLED  LENGTHS  OF  COLUMNS  IN  MRT  TABLE. 

TFNO  hill  CHNTAIN  length  of  longest  column  in  mrt  TARLE. 

DATA  ISTOP/10«20/.  IENO/0/ 

INPUT  CARD  II 

Read  Tme  following  quantities  prom  first  input  caRo— 
nmObNumMER  of  months  POR  WHICH  PlIh  TAPE  IS  TO  BE  MADE 
rnviaBASE  FiELO  OP  VIE* 

nfbnumner  of  multiples  op  POVl  To  re  consiufred 

RMAXaMtXiMIJM  RANGE  IN  KM— CAN  BE  EITHER  10*  OR  ?0* 

rear  Son.  NM0.P0V1.NF.RMAX 
PORmaT(?Ms«PR.O)  I 
IPInF.Ot.S)  NPaS 
no  1 IPal.S 
1 FOVlIPlalPaPOVl 

INPUT  CARO  PI 

RpAO  IN  The  mrt  curve  selector.  tAmoET  SIPt.  target  OIPPERENTIAl 
temperature  (DELTA  T).  ANO  ASPECT  FACTORS I 

ICURvsmMRT  curve  SELECTOR--CHOoSa  1 oF  the  curves  specified  In  data 
statements  in  subroutine  mesOLvE 

SIZEaTARGAT  SIZE— HEIGHT  OP  TOT  (IN  METERS)  SMAU  EST  DIMENSION 
TaRSThTARGET  temp.— delta  OEG,  K oF  IAROET  prom  the  background 
ASPEC'hTO  correct  value  op  ASPECT  RATIO. t. POR  DETECTION  AND 

RECOGNITION.  CORRECT  VALUE  OP  mrt  hY  THE  PACtOR  I/SORTiE/T). 

iNMiiT  hEre  the  value  sorT»e/ti  Pom  oetfction  (Aspect (in  ano 

RECOGNITION  (asPECT(2)),  E IS  Ban  LENOTH-TO-HIOTH  RATIO. 


\\ 

\\ 
I I 

Z'i 

it 

2l 

21 

3< 

«( 

3( 

3I 

41 

41 

41 

Si 


ir<r»nro  ir»riro  irirn  irir*nro  ir.mr  nn  iroiriro  »riro  inr  o 


i 


2M  OET  *SPECT  RATIO,  All  »tc,  RELOA  MRT  IS  CORRECTED 
Bt  dividing  bt  these  input  values, 

REAn  102,  TCURVE,SIZE,TAR6T,ASPtCTa) ,AbPECT(2) 
in?  FORmatiII ,aX,4Flo,AI 

print  lOT,  NMO,FOvltNF,RMA*,ICUHvEtSlTE,T*RSTt*SPPCT 

format (*-INPUT  CARO  VALUES»/*oCarO  I—  NMO«A,I3**  F0V1«»,F5.1  »• 

- nf»«»U,*  RMAXa*,FS,i/*OCAAu  2--  ICUBVE"*, I 3t»  SlZE**tFA, 3,* 

- target  tEhP»»»F6,3,*  AjPCCt  FACtORsi  0£T»*tFT,4,A  REC««,F7,4» 

vritc  hfaoinos  of  Table  to  be  placed  at  thl  beginning  of  output  tape 
WRItF(2,io«)  FoV 

inn  FORMAT! I M| ,AXt*RAN0t«,5XtAAN0LE  IN«tBX,*CVCLES»*2X,5 (3Xt»HRT (•, 

• F4.1 ,•  OEG  FqV)*) ) 

PRlTF(?,l01) 

mi  FORmaT(B*.«IN  KH  mIlLIRAD  per  mm  •,5(6X,»0ETA,7X,«RfC  •)) 

RANGE  CAN  HE  FROM  ,S  TO  10,  KM  IN  STEPS  OF  .5  KM  (RMaX*10,) 

OR  FROM  1.  TO  2O.  KM  IN  STEPS  OF  KM  (RMAAszO.) 

Go  through  Loop  5 oNCE  FoR  each  pAooE 

DO  R la1,20 

01ST«FLOaT(II*RMAX/20, 

FIND  TmF  anGUI.AR  subtense  (THETA)  OF  THE  TARGET  (IN  MILL  IRAQI AN$) , 

AND  Thf  detection  and  hecOQNITIOn  wmITeRIA  (THEIR  SPaTIAL  FRtQS  IN 
ctcles/mrao>  at  a given  range  3v  Calling  surroutine  spatial 

CALI  SPATUL(01ST,S12E*0ET*REC*  THETA) 

CPMiDaOET  S CPM(2)«BEC 
JaO 

GO  THROUGH  LOOP  2 TO  SCALE  MRT  FqR  FOV  OTHtR  THAN  THf  ORIGINAL  MRT 

oaT*  From  a SPECIFIC  GIVEN  SET  OF  IMPACT  mm  I DATA  ANO  CORRECT  IT 

For  aSpaCT  ratio  once  for  E*CH  FqV  (SEE  REPORT  P,  S4«5S) , 

DO  ? IFal,NF 
SCALFaFLOATIlF) 

Then  Once  for  detection  ANO  ONCE  For  pECOGnition 

DO  ? I0Ra1,2 
JaJ*1 

«C»LE  the  MRT  FOR  SOME  FOV  OTHER  Than  THE  ORIGINAL  mrT  DATA 
Py  FALLING  SUBROUTINE  RESOLVE 

CALI  RES«LwE(CPM(IOm) ,1CURVE, scale, MRT, iflaG) 

correct  mrt  for  aspect  ratio  (DIVIDE  HMT  BT  INPUT  CORRECTION  FACTOR 
FnR  aSPFCT-SEF  EXPLANATION  FOR  ASPtCT  qN  INPUT  CARO  p AND  REFER  TO 
RFPORT  P.Sn'SR  ANO  APPENDIX  0), 

RTMp ( J, I 1 bmrt/ASpECT I IOrI 

IFirFLAG,EO,l,AN0,IST0P(J) ,EQ,20)  ISTOPlJIal 

if(iflag.eo,i»ano,istop(J) ,lt,I)  RTMP(j,n«j,g 
IF(TST0P(J> ,GT,1EN0>  lENOalSToPl J> 


•>  CONTiNUe 


16 

c 117 

!;  4PTTE  TABLE  CONTAINING  mPT  0*TA  FOH  EACH  RanGE  AND  EnV  AS  A U8 

t HrAOE«  TO  THE  FIRST  FILE  OF  THE  NEm  FLAP  Tape  (2o  LIncS  TOTAL)  U« 

S:  1?0 

WRlrF(2.inA)  OISTtTHET*tOETt (RTMplKf I) «KbI,J)  1?1 

IftA  FORMAT(6S«F4,l«7XtFS.3t6X«FS*3icitlOFlo*1)  122 

5 continue  123 

C 12* 

i;  KOL  IS  TOTAL  NO,  OF  COLUMNS  OF  PRObABILITItS  TO  BE  bPITTEN  ON  TAPE  125 

5*  126 

KOLaJ  127 

C 1?§ 

L Sn  THROUGH  LOOP  50  ONCE  FOR  EACH  FILE  TO  Rb  WRITTEN  ON  NEW  TAPE  129 

ir  130 

no  FA  MOa),NMO  131 

PEAOli.ijB)  I8LANK  132 

ina  format (//A10)  I33 

C 13A 

!;  FOR  EACH  NEW  FILE  ON  The  NEW  FLIR  TAPE  135 

t WRITE  A HOLtEPlTH  1 ON  next  LINE  OF  THE  new  TAPE  (WHEN  listing  NEW  136 

5.  TAPE  It  skips  to  a new  page— when  heaoInG  the  NFW  TAPF  reads  137 

c A blank  LINA)  130 

!:  13* 

WRITF(2.1 I")  IAO 

lin  FORMATdHi  .10X1  Ul 

C 142 

t RFAO  Transmissions  (previously  calculated  wv  subroutine  lowtran  and  ia3 

t written  On  )APE1>.  144 

U the  transmission  tape  consists  of  unC  file  for  EACH  month  to  be  145 

L processed,  one  months  file  contains  three  lines  of  uEAOInO  or  title  146 

C TmFORMaTIOn  which  are  skipped  over  when  being  read#  followed  by  one  147 

L LlNp  FOR  Each  HOUR  IN  The  MONTH,  UNE  hours  OATA  is  CUMPOSED  of  integer  148 

C values  representing  Month,  day  and  houh,  ano  20  real  values  representing  149 

i The  fractional  transmission  for  RAnGES'OF  .5  TO  10  Km  IN  STEPS  OF  .5  KM  ISO 

L HP  1 To  20  km  IN  steps  OF  1 KM,  131 

L The  format  for  one  HOOHS  data  IS  (2I2,I4»2»*20EA.3» • lS2 

ir  ls3 

10  REAOll.IOOO)  MOH.tRANs  154 

lOAn  FORMATf2T?,I4,2X,20F#..3)  155 

IF(W0F,11  40«1I  ISO 

c 1ST 

t POR  EVERY  HOUPS  data  on  the  new  tame  white  the  FOLLOWING  LINES-.-  1S8 

U • write  * PLANK  line  159 

L a write  P lines  of  hfaOINGS  FOR  c0|  UMNS  OF  THE  TABLE  1 4O 

c a Write  ?n  lines  of  OATA  (ONE  LlNE/RANOt)  161 

i:  162 

11  WRITF(2,111)  MOH  103 

Ml  FORMAT!//, lOXtaMONTH  a •» I2« loXtaQAT  a •, 12, I QX iWuOUR  a WtlA)  164 

wRItE(2,112)  1o5 

112  FORMAT(laX)  166 

WRITE(2,113)  FOV  167 

113  FORMAT(lX,aRANOEa,5X,a  APPARENT  a, lOX.S (F4. 1 ,a  F0va,12X>)  US 

WBtTC(2.114)  169 

114  FORMATdX.aIN  KMa,5A,aTARG£T  TEMPa.S (7X ,aOETa, TX.aHECa) » 170 

^ 

C GO  THROUGH  LOOP  30  ONCE  FOR  EACH  haNgE  (ONE  RANGE  PFH  ROW  IN  TaSLE)  1T2 

ir  173 


ij  4.3  I f k i'  k ..  V— 4-  I 
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00  10  lalflENO 
RANnCaFLUAT ( I ) •RMAX/70* 

4PTMP  1%  apparent  TAROtT  TENPERAtOhE 

APTMPaTPANR ( I ) aTAROT 

8n  THROUGH  LOOP  20  ONCE  FOR  E*CH  COLUMN  OF  pRORARILITIES 

OO  ?0  Jal.KOL 
POErai.n 

IF(RIMP|j,t).EQ.o.o>  00  TO  IS 

NORMALTFEO  rIRNAL  to  noise  ratio  a apparent  TArOET  temp.  (DELTA  T>/MrT 
(CORRECTED  FOP  ASPECT  MV  OIVIOINQ  ov  SQRT(t/7)  ). 

»,  RA  MANNOVEP  REPORT.PART  1 I SiNtE  THE  mht  IS  THAT  VALUE  OF  INCREMENTAL 
TrMPER»TURE  that  PRODUCES  A VaLUc  OF  SiON*L-TO-NOISE  maTIO  SUFFICIENT 
TO  allom  an  orservEr  tu  rreak  Out  the  hrt  test  rar  pattern  at  so  RErcent 
PPORASILITV, (DELTA  T/MhT»  CORMESpOnuS  lO  thF  NORHALIfaTION  EACToR 
RfLATInR  SImnaL-TO-NOISE  to  PRORaSILITT*  OhfRC  delta  T is  the  InCRE- 
mFNtAL  OIFFKaENCE  in  tEmRErAtorE  OtTHFEN  tmF  TArOET  and  Its  OACkOrOuND 

SNR.APTMP/PTHP ( J • I ) 

CALL  SKHROUTINE  CURROB  TO  DETERMINE  THE  RPOHARILITV  oF  DETEC,  (oR  RFCOO.) 
BtVEN  ThC  normalized  signal  to  NOIiE  RATIO 

CALL  CURR0R(8NR»P0ET» 

RROr( JTaMOET 

TR  IF(TRT0R(J» ,LT,I)  PHOR(J)a0.0 

?«  continue 

C 

MR I TF ( 2 , 1 3n ) RANGE  f APTMP • ( PROR ( J» , Ja I ,koL ) 

T3»*  EORMAT(?t,F4.1teXtFS,3.3**10«5*»FS.i) > 

3«  CONTTNUF 
C 

60  TO  10 
AO  ENO  FILE  2 

An  continue 

STOP 

ENO 


irinr  o »cv€rir»€>n  irinrir-rwrir  o 


SU8P0UT1NE  SPATIAL (OISTI t SIZE, OtT»*fC#TMltTAj 


AMBPOUTINC  spatial  FINOS  THE  ANQULaR  SUSTtNSE  (THETa,  in  MILLIRAOS) 
RrOUlREn  FOe  peSOlVINQ  ONE  BAR  IN  ThE  E^^I'^aI-ENT  BAR  CHART  (SEE  REPORT 
P.AR.Sn  FOP  FACH  RANGE  tO  THE  TARpET 

AND  THr  DETECTION  ANO  hfCOBNITION  LhITERIa  (SPATIAL  FHEO,  IN  CYCLES/nRaoj 

convert  Range  TO  heters 
OISt«OIST1*IOOO, 

PTNO  THF  ANGLE  THETA  By  tAKINO  thE  tAN(»ENt  of  HALF  THE  TOt,  SIZE 
(TN  mEtfRSI/RaNOE  (METERS),  DOUBLING  if  To  GET  THE  HhulE  ANGLE, 

ANO  THfn  converting  it  to  mILLIRaOIaNS 

THEta,2. SATAN (SIZE/ (2. POIST) )«l000. 

The  criteria  or  spatial  FREQ.  For  detection  is  1 CYCLE  PER  the 
angle  THETa,  for  recognition  IS  4 CYCLES  FUR  THE  ANGl.F, 

OET«)./THEtA 
RECpn. /THETA 
HFTiipn 
ENO 


B5I 
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iririr-o  iririr  <-»  irinnoo  o i4nrir»r«r'0»rinr  r»  oirirJr»r»ririr»rir»r  or» 


SUBROUTINE  RESOLVE (NESO.ICURVEtSCALirtNRT.IFLtOI 

rurrouttne  rErOlve  inputs  oiven  knumn  NpT  Curves  *no  scale  factors  of  t^e 
spatial  freo.  based  on  FOV.  given  the  spatial  freo.  calculatco  in 
SiiRROUttnE  rPaTIAL*  this  SUBROUTlNt  SCALES  IT  AND  FInUS  THE 
coRRESPONoiNO  MRT  Value  m linear  interpolation  between  two  points  on  the 
otven  mht  curve*  This  approximation  is  ustn  when  the  mrt  is  not 
avatlarie  fhr  a field  of  view  of  Interest  «ut  is  avatlable  for  some 
not  vert  different  sized  field  of  view. 

TCliRVE  INDICATES  mrt  curve  TO  BE  USED 

RfSO  is  the  criteria  for  either  oEi.  oh  RFC.  (CYCLES/mRAO)rSPATiaL 
freouenctrpesolution  of  The  Sensor 

dimension  CURVE(IS*2«S| 
real  mrt 

Input  riven  mot  curve  points  here  in  oata  statements  ••  one  line 

FOR  THE  X coordinates  (SRaTIaL  FHFO.  in  CVCLES/MRAD)  and  one  line 
FOR  THE  T COORDINATES  (MnT  VALUES!  — 2 LiNfcS  FOR  EACH  MRT  CURVE. 

UP  TO  s Mrt  Curves  may  be  entered  *n  the  data  statement  at  ome  time 
HOWEVER  ONLY  ONE  mrt  curve  is  used  PER  RUN.  THE  OESIONATIONS  AT  The 
End  of  each  line  of  the  oata  statement  REFEm  to  which  MRTS  are  available 
In  the  broomam  aI  The  CURRENT  TImE.  TO  SELECT  ONE  OF  THESE  mRT  CURVES 
set  ICUMVE  ROUAL  to  The  position  The  OtSIPtO  MRT  OCCUPIES  IN  THE  LIST. 


data  CURVE/1. tZ. *3. **.tS*.S**7.«HR0.t  T~1  X 

.02. *0ab* *081 1 * II 1 1 .Zpt * 1 * .§*0. * T»1  V 

I * .2* . 3. . A. t S* 1 8* . 7* *8. .Tap . . T*2  x 

.007*.01M*.03<.*0«.Tie2..2..A.1..7PO.*  T-2  y 

I * *2.  .3*  . A * »S*  . 6*  1 7*  .8.  .9.  *8*l*,  . T“3  X 

.002..00A2..008f*Ol3..026..nA2..0A..17s*.S2.6AO..  T-)  Y 

2.  . A*  .6.  .8.  . 1 0.  » 1 0*0.  . H»1  It 

.0082t.0tA..036..0ST..0PM.10Ap..  H.l  y 

1 . 1 1 *B*2* .2* a'2*8. 3*0.3,2. 3* a*3*8*3.8.3*9.3*99.3*0* . PDO  t 

.08.*16t.2..3t.<2..St.A.,7F,.RA«1.08*1.2A.1.32t3A0./  PDO  y 


IFLAOao 

RFSS  Equals  spatial  frequency  at  b«se  fiei  u of  view  acaled  for  the 
multiple  of  FOV  BEING  CALCULATED  At  ThIs  time 
(AFt  p,  ss  mannover  report*  part  U 

REssaREsRarCALE 

look  Up  RESA  along  X-axis  (CTCLES/hhao)— want  to  FINO  corresponoing 
Y Value  (mrt! 

KalruRVF 

IF  resa  (Scaled  spatial  freg.»  is  less  tnan  the  smaliest  x coordinate  (smo 

SET  RESSaSMX  aNO  SET  MMTaSMALLEST  I COURQINaTE  (SHY). 


SMX.CURVFM  »1»K) 
AMVaCUR/En  .2.K) 
IF<aFSS.OF,SMX)  GO  TO  5 
PESsaSHX 


f' 

-5/ 

{ 

, t 
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RCTUMN 


i cnuNT  MOW  many  porwrs  iMe*<e  «we  alums  the  choscn  mrt  cusve. 

L 

M 00  I*  N«l.J5 

tF(C"RVfelM«l*KI .EQ«0.0)  so  TO  11 
in  continue 
NalA 
NNai««l 

C MMa  NO,  SOTNTM  ON  SIVCN  MAT  CORVC  IN  OAT*  WTATEMeNT 
<;  Ms  THE  AOINT  ALONS  THE  SfVEN  MAT  CumvE  rOA  which  AES<  is  SEINO 
i TfSTEn  TO  SEE  ITS  VALUE  LIES  SeTWECN  Tha  VALUE  OF  POINT  N AnO  POINT 
C M.I  ON  THE  H»r  CUAVE 

C 

00  ?n  N»P«NN 

TF(AESS.t*T.CUAVE(N*l«iO)  60  To  Wo 
C 

t at  thin  point  THE  X-COOAOINATE  InTLAVAL  HA5»  rEEn  found,  INTEAPoL*TE  BETwECN 
^ The  Y>COOA01NaTES  To  find  the  exact  MAi  VALUE. 

c 

HOT«CCti»vE(N,2,Ki-CUAVE(N-l,2,KM*(HESS»eUAVE|N-l.l,Kn/ 

• •CUHwE(N,1,K).CuavE(N-1*1,KM*cUAVE(n.I,2,K) 

RETii"N 

?«  continue 
c 

Jr  TF  AESn  CANNHt  SB  foum  along  the  blVEN  MPT  CUAVE  AND  ITS  VALUE  IS  OFF 
c The  CUAVE  Af  THE  HIGH  tNO  THEN  Set  IFLAGal  aNO  SET  HATaTHE  HlwHEST 
k V.COOAoINATA  r>N  THE  GIVEN  MAT  CURVE 
L 

IFLAGal 

HRTaCUAVE (HN«2*K) 

AFTlIHN 

FNO 


rccT  A'/-" 

Dl.j1  1 1 w f ' 
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SUB*«OUT  INC  CUpMOti(SMMtCMPH(il 

c 

C THIS  subroutine  OETERnInES  the  RROyAaiLlTv  OF  DETEC*  (OR  RECUO.) 

C FOR  A UIVEN  StONAL  TO  noISE  RATIO  RV  LOOKTNQ  uR  THE  SNR  VALUE  Un  A 
C cUNDLaTIVE  NORHaL  (OAUSsIAN)  RRObAbILITV  curve  MHERE  the 
C vAluE  UElT*  T/MRTal.O  CORRESPONDS  TO  A RRoRAHILITY  OF  DETECTION  OF 

c SO  Percent 

C (SEE  p*  S3  AND  61  HANNOvER  REPORT*  PaRT  1) 

OImENSION  SRATlOUOl.PROlidO) 

0*Ta  NA/10/ 

INPUT  cumulative  normal  probability  curve  coordinates 

DAT*  S«ATIO/0.**S**6S**M*  1**  1*1*1. .2S*1*S*1*7S*2*/ 

DAT*  PROB/O* **1**2**4**5**6* *75^ *6« *06* 1 */ 

IF(SNR*0T*SRAT10(1>)  go  to  2 
CMPNBaPROBIl) 
return 

2 continue 

NlsNN.1 
DO  3 JaltNi 

if(Onr*oe.sratiou).a*Snr*lT.Sratio(J*1))  Go  To  a 

3 continue 

CMPNBaPROB(NN) 

return 

c 

c interpolate  between  a points  on  N0hM*lI7E0  curve  to  OETeRhINL 
C probability  Op  DET*  (or  REcOO*)  corresponding  To  CAL^CULaTeo 
C SIGNAL  TO  NOISE  RaTiO 

c 

A continue 

XX«SRATIO(K*l>-SRATin(K) 

YYaPROB (K*l I -PROB (K) 

XPa&NR.SRATI0(R) 

CMPHBa ( YY*XP/XX I *PR0R (K> 

return 

ENO 


1C 

2C 

3C 

*C 

SC 

6C 

7C 

n 

lOC 

lie 

12C 

13C 

UC 

19C 

16C 

ITC 

Toe 

20C 

21C 

22C 

23C 

2AC 

29C 

26C 

2TC 

2BC 

aoc 

30C 

31C 

32C 

33C 

3AC 

3Sc 

36C 

3TC 

3BC 

30C 


U 


Lui  i 
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V.  EXAMPLES  OF  OUTPUTS  WHICH  CAN  BE  GENERATED 
FROM  THE  OUTPUT  TAPE  OF  PROGRAM  FLIR 

There  are  many  ways  of  using  the  data  that  has  been  written  on  the  output  tape 
from  Program  FLIR  (hereafter  referred  to  as  the  FLIR  tape).  For  the  most  part 
it  has  been  used  at  IDA  to  generate  graphical  displays  of  the  data.  The  following 
are  some  examples  of  plots  that  can  be  made  with  this  data. 

1.  Prc^abilitv  vs.  Time  Plot 

Figure  2 shows  probability  of  detection  (or  it  can  show  recognition)  plotted  against 
days  of  a given  month.  This  particular  plot  was  done  for  one  FOV  and  two  ranges 
(thus  two  curves).  The  data  necessary  to  make  this  plot  was  extracted  from  the 
hourly  tables  of  the  FLIR  tape.  One  point  per  hour  of  the  month  was  plotted.  The 
hour  of  the  month  was  the  X -coordinate  of  the  point. 

The  Y -coordinate  value  of  the  point  was  determined  by  specifying  the  following 
parameters: 

• File  (or  month)  on  FLIR  tape. 

• FOV. 

• Detection  or  recognition. 

• Range(s). 

The  FOV  and  choice  of  detection  or  recognition  define  which  column  of  the  hourly 
table  (see  Table  2)  the  probability  is  to  be  taken  from.  The  range  specifies  which 
row.  Thus  the  probability  value  occupying  that  position  in  each  hourly  table  for  the 
month  is  extracted. 
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Figure  3 is  a plot  very  similar  to  Figure  2.  In  this  case  transmittcmce  is  plotted 
versus  time.  Because  transmittance  is  plotted  instead  of  probability  it  is  known 
that  the  data  was  extracted  from  the  apparent  target  temperature  column  (trans- 
mission = apparent  target  temperature/  temperature  of  target)  of  the  hourly  table. 
Range  ctgain  was  selected  indicating  the  appropriate  row  of  the  table  (Figure  3 shows 
two  ranges  plotted). 


2.  Ranee  at  Which  Probability  Equals  Specified  Percent  vs.  Time  Plot 


Figure  4 illustrates  another  type  of  data  plotted  against  hours  of  a given  month. 
The  Y-axis  in  Figure  4 represents  the  rcuige  at  which  probability  of  either  detection 
or  recognition  for  a given  FOV  equals  a specified  percent  (e.  g. , 50  percent).  As  in 
Figure  2 the  FOV  and  choice  of  detection  or  recognition  defines  the  appropriate 
column  of  data  in  the  hourly  tables.  However,  the  value  for  range  to  be  plotted  for 
a particular  hour  is  now  determined  by  scanning  the  entire  column  of  probabilities 
and  interpolating  linearly  between  range  values  given  in  the  table  to  get  a range  at 
which  the  probability  is  what  has  beoi  specified  (e.  g. , 50  percent). 


Two  curves  have  been  plotted  in  Figure  4,  one  each  for  two  different  systems. 

To  plot  two  systems  on  one  graph  two  FLIR  tapes  will  have  to  be  read.  These  two 
FLIR  tapes  would  have  been  made  using  the  same  transmission  data  tape  but  inputting 
different  MRT  curves  in  Subroutine  Resolve  of  Program  FUR. 


3.  Probabili 


ified  Percent  at  Given  Ranees  vs.  Time 


Probability  greater  than  or  equal  to  a certain  percent  at  several  ranges  can  be 
plotted  as  a broken  line  plot  to  indicate  hours  of  delay  of  detection  or  recognition 
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in  a given  month.  Figure  5 shows  two  such  plots.  For  the  same  month  for  a given 
FOV  and  five  selected  ranges  one  plot  shows  probability  > 50  percent  euid  the  other 
shows  probability  > 90  percent.  Again  the  FOV  and  choice  of  detection  or  recognition 
indicate  which  column  of  data  to  use  from  the  FLIR  tape  hourly  tables.  Then  for 
each  range  selected  the  probability  in  the  correct  column  is  measured  against  the 
probability  criterion,  for  example  50  percent.  If  the  probability  for  that  range  is 
> 50  percent  a line  is  drawn  for  that  hour,  if  it  is  less  than  50  percent  a blank  space 
is  left  for  that  hour.  The  next  plot  shown  adds  up  the  consecutive  hours  of  delay  and 
presents  the  data  in  a histogram  type  plot. 

4.  Histogram  of  Duration  of  Delay 

Figure  6 is  a histogram  display  of  the  same  type  of  data  presented  in  Figure  5. 
Instead  of  drawing  a line  for  each  hour  where  probability  > specified  percent,  count 
the  number  of  consecutive  hours  where  the  probability  is  less  than  the  specified 
percent  (i.  e. , there  was  a delay).  Then  count  the  number  of  times  a delay  of  that 
duration  occurred  during  the  month.  Again  the  FOV  and  detection  or  recognition 
choices  are  made.  One  histogram  is  plotted  per  range  selected. 

Figure  6 shows  a plot  for  each  of  four  different  ranges.  The  probability  criterion 
in  this  case  was  50  percent. 

5.  Fraction  of  Occurrences  in  Which  Probability  ^ Specified  Percent  vs.  Range 

Figures  7 and  8 show  one  of  the  most  useful  types  of  plots  from  the  FLIR  tape 
data.  Fraction  of  occurrences  of  probability  of  detection  or  recognition  greater  than 
or  equal  to  a given  percent  (for  example,  50  percent)  can  be  interpreted  as  fraction  of 
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Probability  Probability 

>90*  >S<«p 
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successful  events.  Many  sets  of  data  can  be  plotted  versus  range  in  this  manner. 

Figure  7 compares  the  performance  of  different  FLIR  systems  vs.  range  for  a 
particular  month.  To  plot  each  curve  a different  FLIR  tape  was  used  (i.  e. , each  system 
was  on  a separate  tape).  An  FOV  and  detection  or  recognition  are  chosen  to  indicate 
which  column  of  the  hourly  table  contains  the  correct  data.  For  every  range  the 
probability  in  that  column  is  measured  against  the  criterion  probability.  Count  up 
for  every  hour  in  the  month  the  number  of  times  the  probability  for  each  rcmge  is 
greater  than  or  equal  to  the  criterion  probability.  Plot  the  total  number  of  successful 
occurrences  for  each  range  against  the  range. 

Figure  8 shows  basically  the  same  type  plot  as  Figure  7 except  that  it  is  compar- 
ing the  performance  of  one  system  at  two  different  hours  of  the  day  over  a month's 
time.  To  get  the  data  for  this  plot  instead  of  summing  up  the  number  of  times 
probability  > criterion  probability  is  achieved  for  every  hour  of  a month  the  data  is 
summed  once  for  each  6 AM  hour  of  the  month  and  once  for  each  6 PM  hour  of  the 
month. 

Many  different  comparisons  of  the  data  generated  by  Program  FLIR  can  be  made 
using  the  fraction  of  occurrences  versus  range  type  of  plot.  There  are  also  many 
variations  on  the  way  the  FLIR  data  can  be  plotted  versus  time.  The  particular  types 
of  plots  shown  here  are  only  a few  examples  of  the  manner  in  which  the  FLIR  data  can 
be  presented. 
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Figure  7.  FRACTION  OF  SUCCESSFUL  EVENTS  VS.  RANGE  FOR  THREE  SYSTEMS 


O 


Figure  8.  FRACTION  OF  SUCCESSFUL  EVENTS  VS.  RANGE 
FOR  TWO  HOURS  OF  A MONTH 


REFERENCES 


1 


1.  Institute  for  Defense  Analyses,  Effect  of  Weather  at  Hannover.  Federal  Republic 
of  Germany,  on  Performance  of  Electrooptical  Imaging  Systems.  Part  1,  Theory. 
Methodology,  and  Data  Base,  IDA  Paper  P-1123,  L.  M.  Biberman,  August  1976. 

2.  Institute  for  Defense  Analyses,  Effect  of  Weather  at  Hannover.  Federal  Republic 
of  Germany,  on  Performance  of  Electrooptical  Imaging  Systems.  Part  2,  Perform- 
ance of  8.  5 -11  Mm  FLIR  at  Low  Altitude  (U),  IDA  Paper  P-1124,  L.  M.  Biberman, 
March  1977  (CONFIDENTIAL). 

3.  Institute  for  Defense  Analyses,  A Compcurison  of  FLIR  Performance  in  the  3-5  <im 
and  8.  5-11  ^m  Bemds.  Part  3 of  Effect  of  Weather  at  Hannover.  Federal  Republic 
of  Germany,  on  Performance  of  Electrooptical  Imaging  Systems  (^U).  IDA  Paper 
P-1128,  L.  M.  Biberman  and  R.  E.  Roberts,  in  preparation  (SECRET). 

4.  Institute  for  Defense  Analyses,  A Comparison  of  Active  and  Passive  Low-Light- 
Level  Television.  Part  4 of  Effect  of  Weather  at  Hannover.  Federal  Republic  of 
Germany,  on  Performance  of  Electrooptical  Imaging  Systems  (U),  IDA  Paper  P-1202, 
L,  M.  Biberman  and  F.  A,  Rosell,  in  preparation  (CONFIDENTIAL). 

5.  Institute  for  Defense  Analyses,  A Comparison  of  FLIR  Performance  with  8.  5-11  Mm 
Band  with  Active  TV  Performance  in  the  0.  86  Mm  Band.  Part  5 of  Effect  of  Weather 
at  Hannover.  Federal  Republic  of  Germany,  on  Performance  of  Electrooptical 
Imaging  Systems  (U),  IDA  Paper  P-1203,  L.  M.  Biberman  and  M.  L.  Sullivan,  in 
preparation  (SECRET). 

6.  Optical  Physics  Division,  Air  Force  Geophysics  Laboratory,  Atmospheric  Trans- 
mittance From  0.  25  to  28.  5 m:  Supplement  Lowtran  3B  (1976).  AFGL-TR-76-0258, 

J.  E.  A.  Selby,  E.  P.  Shettle  euid  R.  A.  McClatchey,  November  1976. 

7.  F.  A.  Rosell  and  R.  H.  Willson,  "Recent  Psychophysical  Experiments  and  the  Display 
Signal-to-Noise  Ratio  Concept, " Chapter  5 in  L.  M.  Biberman,  ed. , Perception  of 
Displayed  Information.  Plenum  Press,  New  York,  1973,  pp.  167-232. 

8.  Systems  Development  Division,  Defense  and  Electronic  Systems  Center,  Westinghouse 
Electric  Corporation,  Performance  Synthesis  of  Electro -Optical  Sensors.  Report 
AFGL-TR-74-104,  F.  A.  Rosell  and  R.  H.  Willson,  April  1974. 


- 35  - 


